UNCLASSIFIED 


deModuced 
Lif  U» 


ARMED  SERVICES  TECHNICAL  INFORMAnON  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  In  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  Incurs  no  responsibility,  nor  ary 
obligation  \diatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  Implication  or  other¬ 
wise  as  In  any  maimer  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rl^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  In  emy  way  be  related 
thereto. 


i  296  814 

O  - 


<jq 


a.) 
c- 

I '5 

oS5 


■  . 

gmze 

Lj  i  I  ‘  '  '  '  ;  * 

;-i;i 

‘'tS'A 

A.  R«  A.  P, 

AERONAUTICAL  RESEARCH  ASSOCIATES  of  PRINCETON,  INC. 


ARAP  REPORT  NO.  4? 
EXPERIMENTAL  INVESTIGATION  OP  THE 
STRUCTURE  OP  VORTICES  IN  SIMPLE 
CYLINDRICAL  VORTEX  CHAMBERS 


by 

Coleman  duP.  Donaldson 
and 

Richard  S.  Snedeker 


This  research  was  supported  by  the  United  States  Air 
Porce  through  the  Air  Porce  Office  of  Scientific  Research 
of  the  Office  of  Aerospace  Research,  under  Contract  N\imber 
AP  49 (638) -255.  Reproduction  in  whole  or  in  part  is 
permitted  for  any  purpose  of  the  United  States  Government. 


Aeronautical  Research  Associates  of  Princeton,  Inc. 
50  Washington  Road,  Princeton,  New  Jersey 


December  1962 


Abstract 


An  experimental  study  of  the  character  of  vortices  in 
simple  cj'^lindrical  vortex  chambers  completely  open  at  one  end 
is  presented.  Special  attention  is  given  to  the  problem  of 
the  transition  from  one-celled  to  two-celled  vortex  structure. 
It  is  shown  that  for  turbulent  vortices  in  vortex  chambers 
having  length-to -diameter  ratios  in  the  range  0  <  L/D  <  5f 
the  character  of  the  flow  is  primarily  dejjendent  on  the  ratio 
of  a  characteristic  tangential  velocity  to  a  characteristic 
radial  velocity  V/U.  In  addition,  it  is  shown  that  the 
Reynolds  niimber  based  on  the  radial  flow,  while  a  parameter 
of  major  importance  In  laminar  vortices,  plays  but  a  minor 
role  in  turbulent  vortex  motion.  Transition  from  one-celled 
to  two-celled  vortex  structure  started  for  all  vortex 
configurations  tested  (0  <  L/D  <  5)  vjhen  V/U  exceeded 
approximately  3. 
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1 .  Introduction 

Recently,  there  has  been  much  interest  in  vortex  flows. 
This  Interest  has  been  stimulated  by  the  fact  that  there 
have  been  proposed,  within  the  past  fe\i  years  or  so,  a  number 
of  interesting  flo\j  devices  which  depend  for  their  perform¬ 
ance  on  the  properties  of  rapidly  rotating,  axially  symmetric 
flows.  Some  of  the  Interesting  experimental  studies  of 
vortex  motions  that  have  been  made  are  those  of  Kendall 
(1962),  Rosenswieg  (1962),  Ragsdale  (i960),  Donaldson  (I961), 
McCune  and  Williamson  (I961),  Schowalter  and  Johnstone  (i960), 
and  Keyes  (I961). 

It  is  beyond  the  scope  of  this  paper  to  attempt  a 
review  or  critique  of  all  of  these  recent  studies.  Suffice 
it  to  say  that  all  of  them  considered  in  some  detail  the 
effects  of  turbulence  on  vortex  motions  and/or  the  type  of 
£Lxial  flow  patterns  that  exist  in  vortex  chambers  having 
specific  configurations. 

The  present  paper  is  no  exception.  In  this  paper  ’.ve 
shall  deal  with  the  motion  in  simple  cylindrical  vortex 
chambers  which  are  completely  open  at  one  end.  We  shall  be 
interested  not  only  in  the  effects  of  turbulence  on  such 
vortices  but  also  in  the  behavior  and  structure  of  the 
axial  motions  that  exist. 

The  present  work  grew  quite  naturally  out  of  a  rather 
extensive  theoretical  study  of  the  characteristics  of  three- 
dimensional  laminar  vortices  carried  out  by  the  senior  author 
of  this  paper  and  his  colleagues  over  the  past  few  years  (see 
Donaldson  (1956)  and  Donaldson  and  Sullivan  (i960)).  During 
the  course  of  these  Investigations,  theoretical  solutions 
were  found  for  both  one-  and  two -celled  vortex  structures. 
Since  both  one-  and  two -celled  vortex  structures  had  been 
observed  In  the  laboratory  and  there  did  not  exist  In  the 
literature  any  experimental  study  specifically  aimed  at 


1-2 


determining  the  criteria  for  transition  from  one  type  of 
flow  behavior  to  the  other,  the  present  study  was  Initiated 
In  order  to  attempt  to  establish  such  criteria. 

Because  of  the  v;ay  In  which  this  study  originated.  It 
Is  felt  that  It  will  be  helpful  if,  before  describing  the 
experimental  apparatus  used  and  the  experiments  themselves, 
a  very  brief  review  of  the  basic  equations  and  theoretical 
results  of  the  previous  studies  Is  given. 
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2.  Basic  Equations 

In  this  study  ve  will  assxime  that  we  are  dealing  vjith 
an  Incompressible  fluid  of  constant  viscosity.  We  will 
assume  axial  symmetry  so  that  for  a  cylindrical  coordinate 
system  {r,  $,  z)  we  set  the  operator  =  0.  The 

motion  is  further  assumed  to  consist  of  a  steady  part 
with  velocity  components  u,  v,  and  w  and  a  fluctuat¬ 
ing  part  having  components  u  v and  w '  whose 
time  averages  are  zero.  The  time  averaged  equations  of 
motion  for  the  steady  part  of  such  a  flow  satisfy  the 
following  set  of  equations 
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Here,  P  =  p/p,  where  p  Is  the  pressure  and  p  Is  the 
density,  v  is  the  kinematic  viscosity  ii/p,  and  the 
operators  eind  D/D?  are 
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ht 


^  3 
+  u— 

hr 


The  correlations  of  the  form  u'u',  u'v',  etc.  are  the 
negatives  of  the  usual  Reynolds  stress  terms  divided  by  the 
density  (see  Lamb  (1932)  pp.  674-678). 

Some  very  general  conclusions  concerning  this  set  of 
equations  can  be  dravm  from  an  examination  of  their  non- 
dimensional  form.  In  order  to  put  equations  (1)  throui^h  (4) 
in  this  form,  it  is  necessary  to  make  some  assxmiptions 
concerning  the  proper  parameters  by  which  the  Reynolds  stress 
terms  are  to  be  made  non-dimensional.  We  first  note,  however, 
that  in  view  of  equation  (1),  u  and  17  are  closely  connected 
so  that  the  magnitude  of  u  orders  the  magnitude  of  w. 

We  thus  choose  a  reference  radial  velocity  U  and  define 


u  ®  Uu(r,  z) 
w  s  Uw(r,  z) 


(5) 


where  r  =  r/R  and  z  =  z/R  are  non-dimensional  coordinates 
defined  in  terms  of  the  reference  length  R.  Since  the 
magnitude  of  the  tangential  velocity  v  is  not  ordered  by 
either  u  or  w,  v/e  choose  a  second  reference  velocity  V 


for  the  tangential  motion  and  set 
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V  =  Vv(r,  z) 


(S) 


V/e  may  now  return  to  the  question  of  the  non-dimensionali- 
zation  of  the  Reynolds  stress  terms.  V/e  will  assume  that 
these  stresses  are  functions  of  the  related  rates  of  deforma¬ 
tion.  When  account  is  taken  of  the  axial  symmetry,  these 
defonnation  rates  and  the  reference  velocity  on  which  they 
depend  are 
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In  view  of  the  relationships  exhibited  in  equations  (7), 
we  will  make  the  Reynolds  stress  terms  non-dimensional 
according  to  the  following  scheme 
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In  order  to  make  the  pressure  dimensionless,  we  note 
that,  according  to  equation  (2),  a  part  of  the  static 
pressure  must  balance  the  centrifugal  forces  which  are 
proportional  to  the  square  of  the  tangential  velocity.  On 
the  other  hand,  according  to  equation  (4),  a  part  of  the 
pressure  must  account  for  the  acceleration  of  the  fluid  in 
the  axial  direction.  We  will  thus  write 


P  =  V^P^(r,  z)  +  U%(r,  z) 


(9) 


In  terms  of  the  non-dimensional  quantities  defined  above, 
the  basic  equations  (1)  through  (4)  beccmie  (with  t  =  tU/R) 
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(11) 

(12) 

(13) 


We  may  make  the  following  observations  from  the  system 
of  equations  (10)  through  (13).  Firstly,  the  fluid  motions 
are  given  completely  in  terms  of  just  t\io  parameters.  These 
parameters  are  a  Reynolds  nimiber  N  =  UR/v  based  on  the 
reference  radial  velocity  and  a  characteristic  radial 
dimension,  and  the  ratio  U/V,  l.e.  the  ratio  between  the 
characteristic  radial  and  tangential  velocities.  Secondly, 


2-5 


for  high  Reynolds  numbers  when  the  Reynolds  stresses  are 
dominant,  we  would  expect  the  parameter  U/V  to  become 
the  controlling  factor  in  determinin'^  the  motion.  This 
statement  must  be  qviallfled  somev/hat.  V/e  can  never  expect 
the  ratio  U/V  to  become  the  sole  parameter  governing  the 
motion  at  high  Reynolds  numbers  since  the  stress  functions 
a.  •  must  themselves  be  v/eakly  dependent  on  Reynolds  number. 

^  J  p 

In  addition,  the  viscous  term  (v/UR)V  w  in  equation  (13) 
will  turn  out,  in  the  particular  flows  we  will  investigate 
here,  to  behave  in  a  boundary -layer -like  manner  and  hence 
the  Reynolds  number  N  remains  an  important  paraineter 
although,  as  we  shall  see,  it  is  not  the  primary  parameter 
in  governing  the  motions. 

We  may  also  make  the  follov/ing  general  observation 
concerning  the  tangential  velocities  which  are  governed  by 
equation  (12).  We  note  from  this  equation  that  if  both  the 
laminar  and  turbulent  stresses  are  absent  (v  =  0  and 
^ij  “  then  the  product  vr  remains  constant  along  all 
streamlines.  Thus,  if  the  radial  flow  is  Inward,  we  would 
obtain  a  truly  vortex-like  behavior.  That  is,  the  tangential 
velocity  would  Increase  Inversely  v/ith  radius  as  the  flow 
proceeded  towards  the  axis  of  S3rmmetry.  If  only  laminar 
stresses  are  present  (v  0  and  =  0),  the  tangential 
velocity  must  vanish  at  r  =  0  but  the  flow  will  become 
more  vortex-like  the  higher  the  Reynolds  n\4mber  pUR/ia. 

Since  the  mass  flow  into  a  vortex  per  unit  length  in  the 
z  direction  is  proportional  to  pUR,  ;ve  vjould  expect  to 
obtain  more  vortex-like  behavior  of  such  flows  the  higher 
we  made  the  inward  radial  mass  floi;  regardless  of  what  the 
tangential  velocity  is.  Experimentally,  this  type  of 
behavior  is  not,  in  general,  observed.  The  reason  is  that, 
for  the  usual  vorte::  experiments,  the  Reynolds  number  is 
sufficiently  high  so  that  the  vortex  is  turbulent  and  the 
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stresses  which  dominate  the  viscous  behavior  are  the 
Reynolds  stresses.  Turning  our  attention  to  equation  (12) 
for  the  case  of  large  Reynolds  number,  we  see  that  the 
parameter  governing  the  viscous  behavior  of  the  vortex 
Is  U/V.  In  other  v/ords,  for  turbulent  vortices  to  have  a 
more  vortex -like  behavior,  l.e.  tangential  velocities 
approaching  a  l/r  dependence,  v/e  v/ould  want  to  make  the 
ratio  U/V  large,  and  we  would  expect  flows  In  which 
U/V  was  small  to  behave  in  a  manner  somewhat  similar  to 
a  very  viscous  laminar  vortex.  Ue  will  see  in  what 
follows  that  this  rather  general  observation  obtained  from 
dimensional  considerations  is  supported  by  experimental 
results . 
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3.  Boundary  Conditions  and  Existing  Solutions 

In  the  present  experimental  study,  we  are  Interested 
In  the  behavior  of  the  steady  flow  In  a  cylindrical  vortex 
chamber  which  Is  let  Into  an  otherwise  plane  surface  (see 
figure  1).  Let  the  depth  of  this  vortex  chamber  be  L 
and  Its  Inside  radius  R.  The  tani^entlal  motion  Is  Intro¬ 
duced  by  rotation  of  the  porous  cylinder  through  which  a 
radial  flow  Is  maintained  by  virtue  of  a  pressure  differen¬ 
tial  across  the  porous  v/all.  Let  the  tangential  velocity 
at  R  be  V,  where  V  Is  Independent  of  z.  It  Is 
possible  to  arrange  matters  experimentally  so  that  the 
radial  velocity  at  r  =  R  Is  approximately  Independent  of 
z.  Let  this  radial  velocity  at  R  have  magnitude  U. 

At  the  present  time.  It  Is  hopeless  to  seek  a  solution 
of  equations  (10)  through  (13)  that  Is  valid  for  both  the 
region  Inside  and  outside  (above)  the  vortex  cup  Just 
described.  The  best  that  can  be  done  Is  to  seek  solutions 
to  related  problems  having  far  simpler  boundary  conditions 
and  to  use  these  solutions  to  drav;  Inferences  concerning 
the  behavior  of  the  more  general  problem. 

Since,  experimentally,  one  finds  that  what  happens  In 
the  region  exterior  to  the  vortex  cup  depends  on  the  type 
of  flow  produced  within  the  vortex  cup.  It  would  seem 
logical  to  seek  solutions  that  would  be  approximately  valid 
within  the  vortex  chamber.  Even  this  simplified  problem  Is, 
In  general,  very  difficult  because  of  the  behavior  of  the 
viscous  boundary  layers  on  the  floor  of  the  vortex  cup.  If 
these  effects  are  neglected  we  obtain  a  problem  which,  at 
least  for  the  case  of  launinar  flow,  is  tractable  from  the 
point  of  view  of  analysis. 

Following  this  line  of  reasoning,  if  we  neglect  the 
end  effects  of  both  the  open  and  closed  ends  of  the  vortex 
chamber,  we  should  seek  a  steady  solution  to  equations  (10) 


Figure  1.  Vortex  chamber  geometry. 
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through  (13)  subject  to  the  follovdng  boundary  conditions 


u(l,  z)  =  1 
v(l,  z)  =  1 
w(l,  z)  =  0 
w(r,  0)  =  0 


u(0)  =  0 
v(0)  =  0 

5w\  =  0 

'dr/r=0 


(14) 


For  the  case  of  laminar  flow  (v  0,  a..  ^  =0),  such  motions 
have  been  studied  extensively  by  Donaldson  (1955)  and  by 
Donaldson  and  Sullivan  (190O).  The  motions  are  svcn  that 
the  radial  and  tangential  velocities  are  independent  of  z 
while  the  axial  velocity  is  of  the  form  w  =  zf(r).  The 
pressures  and  P2  are  of  the  form 


Pi  =  Pi(r) 
Pg  =  Cz^ 


(15) 


where  C  is  a  constant.  Under  these  conditions,  equations 
(10)  and  (13)  become  independent  of  (11)  and  (12)j  that 
is,  the  radial  and  axial  motions  become  independent  of  the 
tangential  motion.  Once  equations  (10)  and  (13)  are 
solved  for  u(r)  and  f(r),  equation  (12)  may  be  solved  for 
the  tangential  motion,  v  =  v(r).  Finally,  equation  (11), 
from  which  the  term  containing  P^  disappears,  may  be 
solved  for  the  pressure  P^  which  is  associated  \^rith  the 
balancing  of  certrifugal  forces,  radial  accelerations,  and 
stresses  acting  in  the  radial  direction. 

The  most  interesting  single  feature  of  the  results  of 
Donaldson  and  Sullivan's  study  which  is  of  importance  here 
is  the  fact  that  they  found  that  for  inward  directed  radial 
flows  there  existed  two  distinct  solutions  for  the  motion 
in  the  vortex  chamber  corresponding  to  either  a  one-  or 


3-3 


two -celled  flow  configuration  as  illustrated  in  figure  2. 
There  is  nothing  in  the  solutions  themselves  to  Indicate 
which  of  these  tv/o  motions  is  preferred.  However,  given  a 
vortex  chamber  of  finite  length  L  such  as  is  shown  in 
figure  1,  one  can  reason  what  might  determine  the  transi¬ 
tion  from  one-  to  tvjo-celled  flov/. 

Consider  the  vortex  chamber  and  the  flow  emanating 
from  it  shown  in  figure  3.  For  zero  or  very  small  tangen¬ 
tial  velocities  of  the  chamber  v;all,  the  streamlines  are 
very  much  as  shown.  The  flow  leaves  the  vortex  chamber 
almost  parallel  to  the  axis  of  s^nranetry  and,  as  it  does  so, 
it  entrains  by  viscous  action  a  certain  amount  of  the  gas 
in  the  region  into  v;hich  it  flows.  The  pressure  Just  out¬ 
side  the  exit  of  the  vortex  chamber  at  point  A  will  be 
equal  to  ambient  pressure  far  from  the  vortex  chamber.  At 
the  point  B  on  the  centerline  in  the  exit  plane,  the  static 
pressure  will  differ  from  that  at  A  by  an  amovint  which  can 
be  found  by  integrating  equation  (11)  vflth  respect  to  r 
from  the  centerline  to  the  full  radius  R.  For  zero  rotation 
the  static  pressure  at  B  is  slightly  higher  than  the  pres¬ 
sure  at  A  due  to  the  decelleration  of  the  radial  flov;. 

As  the  chamber  wall  is  rotated  faster  and  faster,  the 
pressure  at  B  drops  rapidly  due  to  the  necessity  of 
balancing  the  centrifugal  forces.  Obviously,  for  a  given 
V,  this  pressure  drop  will  be  larger  the  more  vortex-like 
the  tangential  velocity  profile  becomes.  Now  let  us  consider 
the  total  pressure  at  B.  The  total  pressure  at  B  will  be 
made  up  of  the  static  pressure  at  B  plus  the  dynamic  pres¬ 
sure  due  to  the  axial  component  of  velocity  along  the  center- 
line.  Since  this  axial  flow  at  the  exit  plane  of  the  vortex 
chamber  must  account  for  all  the  flo\;  that  entered  through 
the  sides  of  the  chamber,  the  average  velocity  in  the  exit 
plane  w  must  be  given  by 

irR^w  =  2TrRLU 

w  =  2  A 
R 


or 


On* -celled  flows  Two-celled  flows 


Figure  2.  Typical  tangential  and  axial  velocities  In  one- 
and  two-celled  vortices  after  Donaldson  and  Sullivan  (i960). 


Figure  3.  Schematic  view  of  streamline  pattern  for 
nonrotating  vortex  chamber. 
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Thus  for  a  given  radial  flow  U  v;e  would  expect  the 
dynamic  pressure  at  the  point  B  to  be  proportional  to 
the  square  of  the  length-to -diameter  ratio  (L/D)  of  the 
vortex  chamber. 

Now  for  a  given  chamber,  let  us  follow  what  happens 
to  the  total  pressure  at  the  point  B  as  the  vortex 
chamber  is  taken  from  zero  to  a  high  rotational  velocity. 
Initially,  with  V  =  0,  the  total  pressure  is  higher  than 
ambient  pressure  by  an  Increment  that  is  made  up  of  t;7o 
parts;  a  static  part  arising  from  the  terms  containing  u 
in  equation  (10)  and  a  dynamic  part  arising  from  the  mass 
flow  leaving  the  chamber.  VJhen  the  vortex  chamber  is 
rotated,  the  static  pressure  at  B  drops.  It  may  in 
fact,  for  a  given  depth  of  chamber,  drop  so  much  that  the 
total  pressure  at  B  is  less  than  the  ambient  pressure 
into  which  the  vortex  chamber  is  exhausting.  In  such  a 
case,  the  one -celled  flow  we  have  been  envisioning  cannot 
exist.  We  must  suppose  Instead  that  a  two-celled  flovj 
exists.  Actually,  due  to  the  finite  length  of  a  real 
vortex  chamber  and  the  attendant  end  effects,  this  transi¬ 
tion  from  one-  to  two-celled  flow  does  not  occur  suddenly. 
As  we  shall  see  from  the  experimental  results  presented 
later,  as  the  velocity  of  rotation  is  Increased,  a  second 
cell  starts  well  above  the  vortex  chamber  and  is  gradually 
sucked  into  the  center  of  the  vortex  chamber  Itself. 

These  ideas  may  perhaps  be  somewhat  clearer  if  we 
present  a  few  numerical  results  based  on  the  theoretical 
solutions  of  Donaldson  and  Sullivan.  In  figure  4  there  is 
plotted  as  a  fvinction  of  V/U  the  non-dimensional  pressure 
difference  between  the  static  pressures  at  the  exit  plane 
and  on  the  centerline  (point  B)  of  a  vortex  chamber  having 
a  length-to -diameter  ratio  of  five.  As  mentioned  before, 
two  parameters  (V/U  and  N)  govern  the  motion  so  that  the 
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Reynolds  number  N  appears  as  an  Independent  parameter. 

As  expected  the  static  pressure  drops  off  as  the  square  of 
V/U  In  order  to  balance  the  centrifugal  forces.  Also 
as  expected  for  each  family  of  vortices  the  pressure  drop 
Is  more  rapid  the  higher  the  numerical  value  of  the 
Reynolds  number  (the  Resmolds  number  N  is  negative  by 
convention  since  Inv/ard  velocities  carry  a  negative  sign). 
Also  Included  on  figure  4  for  comparison  is  the  curve 
showing  the  static  pressure  drop  that  would  be  experienced 
If  the  vortex  were  not  to  turn  In  a  vortex-llke  manner  at 
all,  but  were  to  have  infinite  viscosity  and  hence  turn 
within  the  vortex  chamber  as  a  solid  body. 

The  critical  value  of  V/U  for  which  the  total 
pressure  at  point  B  of  a  vortex  chamber  is  Just  equal  to 
the  ambient  pressure  is  plotted  in  figure  5  as  a  function 
of  the  length -to -diameter  ratio  of  the  chamber.  Here 
again  the  Reynolds  number  N  is  an  independent  parameter. 

As  one  might  expect,  the  lower  the  Reynolds  number  the 
higher  must  be  the  V/U  of  the  vortex  chamber  to  achieve 
the  critical  pressure  drop.  A  very  interesting  fact  in 
regard  to  these  results  is  the  extremely  low  ratio  of 
V/U  that  is  required  to  make  a  one -celled  flow  an 
impossibility  for  normal  vortex  chambers. 

This  is  about  all  the  information  concerning  the 
behavior  of  laminar  solutions  that  it  will  be  useful  to 
present  at  this  time.  Since  we  v;ill  be  primarily  interested 
in  turbulent  vortices,  we  will  proceed  to  a  short  discussion 
of  how  we  might  expect  the  results  discussed  above  to  be 
modified  for  the  case  of  turbulent  motion. 
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4.  Turbulent  Vortices 

At  the  present  time,  far  too  little  Is  known  concern¬ 
ing  the  nature  of  turbulent  stresses  in  general  or,  for 
that  matter,  about  turbulent  stresses  in  the  particular 
case  of  the  vortex  chambers  under  discussion  here  to  permit 
an  analytical  solution  to  equations  (lO)  through  (13) • 
However,  if  one  ass\.imes  a  typical  mixing  length  formulation 
for  the  turbulent  stresses  either  of  the  form 


and  similar  expressions  for  the  other  stresses,  it  is  not 
hard  to  show  that  the  very  simple  flows  of  the  form 

u  =  u(r) 

V  =  v(r)  (13) 

w  =  zf(r) 

obtainable  in  the  laminar  case  can  no  longer  identically 
satisfy  the  equations  of  motion.  Because  of  this  difficulty, 
as  well  as  the  conceptual  difficulties  inherent  in  applying 
a  mixing  length  or  other  empirical  approach,  no  theoretical 
analyses  similar  to  that  given  by  Donaldson  (I961)  and 
pursued  by  McCune  and  V/llllaimson  (I96I)  have  been  attempted. 

Actually,  in  practice,  the  turbulent  motions  observed 
can  come  very  close  to  satisfying  the  conditions  given  in 
equations  (18).  Thus,  from  a  practical  point  of  view. 
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some  insight  into  the  behavior  of  turbulent  vortices  can 
be  gotten  by  considering  them  to  behave  in  a  manner  similar 
to  laminar  vortices.  In  this  case,  however,  to  describe 
the  viscous  effects  we  replace  the  Reynolds  nvimber  UR/v 
by  a  number  proportional  to  U/V.  That  is  for  turbulent 
vortices,  we  take  an  effective  Rej'-nolds  number  N  defined 
as 


N*  =  k  -  (19) 

V 

As  mentioned  before,  this  reduces  the  number  of 
parameters  governing  the  vortex  motions  under  consideration 
to  Just  one,  namely,  U/V.  What,  then,  would  we  expect 
to  find  in  regard  to  the  behavior  of  such  parameters  as 
are  plotted  in  figures  4  and  5  for  the  case  of  turbulent 
vortices.  Since  there  is  only  one  parameter,  the  families 
of  curves  dependent  on  Reynolds  number  for  both  the  one- 
celled  and  two-celled  flows  plotted  in  figure  4  would 
collapse  into  Just  one  curve  for  each  family.  Further, 
since  the  critical  rotational  speed  (V/U)  as  shown  in 
figure  5  is  no  longer  a  function  of  Reynolds  number,  we 
would  expect  to  find  that  (V/U)  is  a  function  of  L/D 

w  X 

alone.  In  view  of  this  fact,  we  v/ould  expect,  for  a  given 
vortex  chamber,  that  a  pressure  plot  such  as  that  shown  in 
figure  4  would  also  consist  of  Just  a  single  curve.  This 
curve  would  start  out  from  V/U  =  0  following  the  single 
curve  of  pressure  versus  V/U  for  one-celled  vortices.  As 
it  approached  some  critical  value  of  V/U,  the  curve  would 
transit  so  as  to  follow  the  curve  of  pressure  versus  V/U 
for  two-celled  vortices.  While  following  the  pressure 
curve  for  any  one  particular  type  of  flow  (one-celled  or 
two-celled),  we  would  expect  these  curves  to  have  less 
negative  slope  than  the  curves  of  constant  Reynolds  number 
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since  we  are  effectively  Increasinj^  the  viscous  action  or 
lowering  the  effective  Reynolds  number  as  V/U  is  increased. 
Before  passing  on  to  a  discussion  of  experimental 
results,  it  is  worthwhile  to  point  out  an  interesting 
constraint  on  the  total  physical  stress  in  a  vortex, 

whether  it  is  laminar  or  turbulent,  for  the  special  case 
for  which  the  tangential  velocity  is  a  function  of  the 
radial  coordinate  alone.  In  this  case,  the  equation 
governing  the  tangential  velocity,  equation  (12),  may  be 
written  in  terms  of  the  total  stress 

as 

(21) 

dr  ^  dr 


Integrating  equation  (21)  and  noting  that  =  0  at 

f  =  0  one  obtains 


(ru)dr  (22) 

or 

Thus  the  total  stress  and  hence  the  Rejmolds  stress  can  be 
determined  from  a  measurement  of  the  mean  velocity  profile  s 
u  and  V.  This  relationship  is  often  useful  in  considering 
the  flow  within  turbulent  vortices  (see  for  example  Donaldson 
(1961)  and  McCune  and  Uilllamson  (I96I)). 
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5.  Experimental  Apparatus 

In  order  to  check  experimentally  the  nature  of  the 
cellular  structure  of  simple  vortices  such  as  have  been 
discussed  above,  an  apparatus  was  constructed  that  was 
designed  specifically  to  be  able  to  test  the  conjectures 
put  forward  In  the  preceding  sections. 

5.1.  Vortex  chamber 

The  chamber  used  to  produce  the  desired  three-dimension¬ 
al  vortex  flows  was  designed  so  as  to  provide  the  greatest 
flexibility  In  establishing  the  various  boundary  conditions 
of  Interest.  The  vortices  were  produced  Inside  a  rotating 
tube  of  uniformly  porous  ceramic  material  11-1/2  Inches 
long,  with  an  Inside  diameter  D  of  2.2  Inches  and  a 
3/l6-lnch  wall  thickness.  This  tube  was  supported  vertical¬ 
ly  between  end  plates,  which.  In  turn,  were  mounted  on  the 
Inner  races  of  a  pair  of  ball  bearings .  The  bearing  outer 
races  were  held  In  fixed  end  plates  which  also  acted  as 
end  covers  for  the  steel  pleniam  cl:amber  surrounding  the 
rotating  assembly  (see  figure  6).  A  highly  effective 
labyrinth  seal  was  machined  In  each  rotating  and  stationary 
end  plate  to  minimize  the  leakage.  The  tube  was  rotated 
by  means  of  a  belt -driven  pulley  attached  to  the  lower 
rotating  end  plate.  A  3/^-horsepov;er  direct  current  motor 
was  used  to  drive  the  assembly.  The  power  supply  was  a 
variable  voltage,  rectifier  type  vAilch  made  possible 
rotational  speeds  up  to  2400  rpm.  For  the  2.2-lnch  tube 
used,  a  maximum  tangential  wall  velocity  of  about  24  ft/sec 
could  thus  be  obtained.  Air  was  supplied  from  a  100-cublc 
foot  storage  tank  at  200  pslg.  An  automatic,  spring-loaded 
regulator  valve  v^as  xised  to  maintain  the  desired  stagnation 
pressure  In  the  plenuir.  chamber.  The  highest  working  pressure 


Rotating 

chamber 
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generally  used  (about  3  psig)  was  found  to  give  a  radial 
velocity  at  the  wall  of  about  6,5  ft/sec  which  corresponds 
to  a  radial  Reynolds  number  N  =  UR/v  «  -3800. 

The  effective  length  L  of  the  vortex  "cup"  was 
regulated  by  means  of  a  movable  floor  plate  Inserted  at 
the  lower  end  of  the  tube.  The  axial  position  of  this  plate 
was  adjusted  by  means  of  a  motor-driven  Jack  screw.  The 
diameter  of  the  plate  was  slightly  less  than  the  inside 
diameter  of  the  tube  so  that  no  contact  with  the  rotating 
tube  was  possible.  The  slight  gap  at  the  edge  of  the 
plate  was  Judged  to  have  little  effect  on  flow  conditions 
in  the  vortex  cup.  The  maximum  ratio  of  length  to 
diameter  L/D,  using  the  2.2-lnch  tube,  was  slightly  over 
5. 

The  entire  apparatus --vortex  chamber,  regulator,  drive 
motor,  and  Jack  screw  assembly — v;as  moiinted  in  a  rigid 
steel  frame  which  provided  convenient  mounting  locations 
for  necessary  instrumentation  (see  figures  7  and  8). 

The  device  Just  described  was  thus  capable  of  estab¬ 
lishing  vortex -type  flows  within  the  porous  tube  for  any 
tangential  velocity  up  to  24  ft/sec,  any  radial  Reynolds 
nwiber  up  to  -3800,  and  for  any  length-to-diameter  ratio 
up  to  5j  with  all  three  parameters  Independently  variable. 

5.2.  Instrumentation 

The  radial  pressure  distribution  at  the  bottom  of  the 
vortex  was  measured  by  means  of  six  static  pressure  holes 
in  the  floor  plate.  Static  pressure  and  total  pressure  on 
the  center  line  for  various  axial  locations  (0  <  z/L  =  1.0) 
were  measured  by  means  of  suitable  tubes  made  of  .032 -inch 
O.D.  stainless  steel  tubing  and  mounted  either  at  the  center 
pressure  tap  of  the  floor  plate  or  from  above,  ''?I:e  most 
reliable  performance  v/as  achieved  vjlth  a  dual  probe  having 


Figure  7*  Vortex  chamber  mounted  In  stand  showing  drive 
motor  and  chamber  floor  actuating  assembly.  The  variable 
voltage  direct  current  power  supply  Is  at  the  right. 


Figure  8.  View  of  vortex  apparatus  showing  porous  cylinder 
and  5-hole  prism  type  probe  In  position. 
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separate  static  and  total  pressure  tips  each  of  which  could 
be  positioned  in  the  region  of  Interest  by  rotating  the 
tube  90°  about  its  axis — the  probe  axis  being  normal  to 
the  vortex  axis.  This  probe  was  mounted  above  the  chamber 
and  its  use  was  limited  to  positions  near  the  top  of  the 
vortex  cup. 

Velocity  profile  measurements  were  made  with  two 
commercially  available,  pitch  and  yaw  sensitive,  calibrated 
probes.  Each  was  of  the  5-hole  type,  one  having  a  prism 
tip  and  the  other  a  conical  tip.  The  diameter  of  each 
probe  at  the  tip  was  .12 -inch.  (Factors  influencing  the 
use  of  such  probes  in  shear  flows  are  discussed  in 
Section  6.)  The  prism  probe  was  used  for  velocity  distribu¬ 
tions  both  at  the  top  of  (z/L  =  l.O)  as  well  as  Inside  the 
vortex  cup.  The  conical  probe  was  used  only  at  the  top. 

All  pressures  \7ere  read  from  an  adjustable -tilt 
alcohol  manometer  consisting  of  four  U-tubes  suitably 
connected  to  the  probes. 


5.3.  Calibrations 

Measurement  of  the  supply  tank  blowdown  rate,  assuming 
constant  tank  temperature,  was  used  to  determine  the  approxi¬ 
mate  mass  flow  rate  as  a  fvinctlon  of  plentun  chamber  pressure. 
The  leakage  flow  through  the  labyrinth  seals  was  determined 
with  a  solid  tube  in  place  of  the  porous  tube,  and  this  flow 
was  subtracted  from  the  over-all  mass  flow  to  determine  the 
net  flow  through  the  porous  wall.  The  leakage  rate  was 
found  to  be  less  than  1  per  cent  of  the  total  for  all 
conditions  tested.  The  over-all  mass  flow  was  Independent 
of  both  rotational  speed  and  length-to -diameter  ratio  of 
the  vortex  cup,  so  that  a  single  calibration  curve  was 
sufficient  for  all  possible  conditions.  Values  of  U  and 
N  were  based  on  this  curve  assvimlng  a  uniform  axial  distribu¬ 
tion  of  mass  flow  through  the  porous  wall.  Rotational  speed 
was  calibrated  as  a  function  of  drive  motor  armature  voltage 
by  means  of  a  Strobotac. 
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6,  Experimental  Results 

Since  the  principal  objective  of  this  work  was  the 
establishment  of  criteria  for  determining  the  type  of 
cellular  structure  of  vortices,  procedures  for  detecting 
the  existence  of  such  structure  were  first  developed. 
Assximlng  that  for  a  radial  Inflow  the  only  possible 
structures  have  either  one  or  two  cells,  as  Indicated  by 
the  laminar  theory,  two  methods  of  detecting  such  structure 
were  tried;  (1)  measurement  of  static  and  total  pressures 
on  the  center  line  and  (2)  direct  measurement  of  component 
velocity  profiles. 

6.1.  Static  and  total  pressures  on  the  center  line 

As  defined  herein,  a  two  cell  vortex  with  radial  Infloiv 
at  the  outer  boundary  (porous  wall)  must  consist  of  a 
peripheral  region  of  upward  axial  flow  (directed  out  of  the 
top  of  the  vortex  cup)  and  a  central  region  of  downv;ard 
axial  flow  (directed  Into  the  vortex  cup).  The  upward 
axial  flow  region  thus  consists  of  fluid  that  entered 
radially  at  the  outer  boundary  plus  the  fluid  that  entered 
axially  at  the  center  and  ultimately  reversed  Its  direction. 
(For  an  axially  symmetric  flov;,  of  course,  radial  and 
tangential  velocity  components  must  vanish  at  the  center.) 
Detection  of  a  two -celled  flow  was  thus  based  on  the 
presence  of  either  a  zero  or  negative  (downward)  axial 
velocity  on  the  center  line — zero  axial  velocity  Indicating 
the  lower  extremity  of  the  center  cell. 

Total  head  tubes  movinted  at  the  bottom  of  the  vortex 
cheimber  facing  upward  In  an  axial  direction  and  tubes 
mounted  from  above  facing  downward  were  used  In  conjunction 
with  similarly  mounted  static  tubes  to  make  extensive 
pressure  surveys  for  a  range  of  radial  Reynolds  niimber  II, 
tangential  wall  velocity  V,  length-to -diameter  ratio  L/D, 
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and  centerline  locations  z/L.  Equality  of  total  and 
static  pressures  so  measured  under  Identical  boundary 
conditions  was  assumed  to  be  indicative  of  zero  axial 
velocity.  Consistent  results,  however,  were  difficult  to 
achieve  in  this  v\ray  because  the  necessary  exact  matching 
of  boundary  conditions  for  corresponding  total  and  static 
pressure  measurements  made  during  separate  iTuns  was  extreme¬ 
ly  difficult.  Subsequent  use  of  the  dual  probe  described 
earlier  resulted  in  more  reliable  data  since  the  total  and 
static  pressures  at  the  same  point  were  measured  at  very 
nearly  the  same  Instant.  This  probe  was  used  to  measure  the 
outward  axial  component  of  total  pressure  as  well  as  the 
static  pressure  on  the  center  line  at  z/L  =  1.0  for 
various  values  of  the  parameters. 

With  L/D  and  N  held  constant,  readings  were  taken 
for  a  niimber  of  tangential  velocities  V,  Typical  results 
are  shown  in  figures  9,  10,  11,  12,  and  13  for  several 
radial  Reynolds  nianbers  and  several  depth-to -diameter 
ratios.  The  data  are  presented  in  the  form  of  a  pressure 
coefficient  based  on  radial  dynamic  pressure  at  the  v/all 
plotted  against  the  non-dimensional  velocity  ratio  parameter 
V/U.  For  very  small  tangential  velocities,  it  is  seen 
that  the  difference  between  the  two  curves  represents  the 
dynamic  pressure  due  to  the  upward  (positive)  axial  velocity 
on  the  center  line.  As  V  is  Increased,  there  is  a  drop 
not  only  in  the  static  and  total  pressures,  but  also  in 
the  difference  between  them.  This  convergence  is  Indicative 
of  a  decreasing  upward  axial  velocity  on  the  center  line. 
Finally,  the  two  pressures  become  equal,  at  which  point 
the  upward  axial  velocity  is  apparently  zero.  Further 
increase  of  V/U  beyond  (V/U)  shows  evidence  of  a 

downward  axial  velocity.  With  the  total  head  tip  of  the 
probe  still  facing  downward,  it  begins  to  register  pressures 
slightly  but  consistently  less  than  does  the  static  tip. 


Figure  9.  Total  and  static  pressures  at  the  center  of  the 
exit  plane  (point  B)  as  functions  of  V/U. 


Figure  10.  Total  and  static  pressures  at  the  center  of  the 
exit  plane  (point  B)  as  functions  of  V/U. 


Figure  11,  Total  and  static  pressures  at  the  center  of  the 
exit  plane  (point  B)  as  functions  of  V/U. 


This  is  what  one  mi^ht  expect  if  the  total  head  tube  v/ere 
now  considered  to  be  sensing  a  base  pressure  due  to  the 
axial  flow  directed  away  from  its  open  end.  This  condition 
persists  to  the  highest  V/U  possible.  The  few  cases 
where  individual  points  are  reversed  all  occur  in  the 
vicinity  of  the  peculiar  region  of  sudden  pressure  rise 
and  are  therefore  attributed  to  the  difficulty  of  maintain¬ 
ing  steady  flow  conditions  in  this  region. 

In  all  the  above  cases,  for  a  given  cup  length-to- 
dlameter  ratio  L/D,  the  value  of  (V/u)  was  found  to  be 
virtually  Independent  of  radial  Reynolds  number  over  the 
range  tested.  The  dependence  of  (V/U)  on  L/D  Is  shown 
in  figure  l4  with  the  values  of  N  indicated.  The 
significance  of  the  various  features  of  figures  9  through 
14  is  discussed  in  section  7. 

6,2.  Velocity  profiles 

Typical  tangential  and  axial  velocity  profiles  for 
one  of  the  cases  just  discussed  (N  =  -900,  L/D  =  5)  are 
shown  in  figure  15.  These  profiles  v;ere  determined  by 
using  the  5-hole  conical  probe  mounted  above  the  vortex 
chamber.  Certain  considerations  involving  probe  effects 
on  the  accuracy  of  these  profiles  are  discussed  later  in 
this  section. 

As  v/u  is  increased,  a  decrease  in  outward  axial 
velocity  near  the  center  is  immediately  apparent  while  the 
outer  portions  of  the  profile  are  essentially  unchanged. 

The  tangential  profile  at  first  exhibits  increasing  veloci¬ 
ties  near  the  center,  but  the  maximum  value  soon  begins  to 
drop  with  increasing  V/U  while  at  the  same  time  the  radial 
location  of  the  maximum  is  displaced  outward.  These  trends 
continue  until  an  inv/ard  axial  velocity  on  the  center  line 
approximately  equal  to  the  maximum  outward  value  is 
achieved.  For  higher  V/U,  the  maximum  Inward  velocity 
changes  only  slightly  while  the  radial  extent  of  the  inner 
cell  increases  until  it  finally  occupies  a  large  portion 


Figure  14.  Variation  with  L/D  of  measured  values  of 
(V/U)cr  point  B  (z/L  =  1). 


Figure  1^.  Tangential  and  axial  velocity  profiles  measured 
in  the  exit  plane  with  the  conical  tip  probe,  Inviscid  and 
solid-body  rotation  tangential  profiles  are  shown  by  heavy 
dashed  lines. 
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of  the  2-celled  system. 

Since  it  was  not  possible  to  measure  negative  w 
directly  with  the  5-bole  conical  probe,  those  values  shown 
on  the  center  line  are  computed  from  the  static  pressures 
measured  with  the  dual  probe,  assuming  a  total  pressure 
equal  to  atmospheric.  Positive  values  of  w  on  the  center 
line  measured  with  both  probes  shov;ed  excellent  agreement 
except  in  the  sensitive  region  near  vv  =  0. 

6.3.  Probe  effects  on  velocity  measurements 

By  far  the  most  significant  source  of  experimental 
error  is  thought  to  lie  in  the  use  of  the  5-hole  probes 
in  regions  of  high  shear  and  near  solid  boundaries. 

A  description  of  the  operation  of  the  5-hole  conical 
probe  is  now  given  in  order  to  clarify  the  discussion  of 
such  probe  error  sources.  This  probe  had  a  single  "total 
head"  hole  at  the  apex  of  the  conclal  tip  with  4  "static" 
holes  a  short  distance  down  the  conical  surface,  the 
latter  being  spaced  90°  apart  and  equidistant  from  the  apex. 
One  opposed  pair  of  static  holes  was  aligned  laterally  v;ith 
respect  to  the  probe  shaft  and  the  other  pair  longitudinal¬ 
ly.  The  probe  was  mounted  with  respect  to  the  vortex 
chamber  as  shown  in  figure  I6.  It  could  be  traversed 
radially  in  either  of  two  directions  so  that  either  pair 
of  holes  could  be  in  radial  or  tangential  alignment.  In 
order  to  determine  the  local  flov/  direction,  at  a  particular 
radial  location,  the  probe  was  rotated  about  its  shaft  until 
equal  pressures  were  registered  by  the  lateral  holes.  In 
this  position,  the  probe  lay  in  a  plane  containing  the 
apparent  velocity  vector.  The  angle  between  the  probe  and 
the  velocity  vector  was  then  determined  from  the  calibrations 


Figure  l6.  Schematic  views  of  conical  tip  probe  mounting 
arrangement.  Transverse  travel  (traverse)  Is  Indicated  on 
the  left  and  axial  travel  on  the  right. 


by  measuring  the  pressure  difference  between  the  longitudinal 
holes.  Calibrations  also  gave  the  true  total  pressure 
corresponding  to  the  "total  head"  hole  reading  for  that 
angle.  True  dynamic  pressure  was  similarly  calibrated  for 
the  difference  betv/een  lateral  hole  "static"  pressure 
readings  and  "total  head"  hole  readings  as  a  function  of 
the  angle.  V/lth  the  magnitude  and  direction  of  the  velocity 
vector  so  determined,  the  components  were  computed.  vJhen 
used  in  this  way  in  flows  with  regions  of  high  shear,  it 
is  clear  that  misleading  data  may  easily  result  since 
holes  on  opposite  sides  of  a  symmetrical  body  Immersed  in 
a  shearing  flow  will  show  a  pressure  difference.  Near 
solid  boundaries,  too,  there  is  a  distortion  in  the  floi; 
pattern  around  the  probe  due  to  an  Interference -type  effect 
which  also  changes  the  angle  at  which  the  flovj  approaches 
the  probe.  A  third  consideration  is  the  radial  static 
pressure  gradient  balancing  the  centrifugal  force.  The 
finite  radial  distance  between  the  radially  oriented  pair 
of  holes  may  be  sufficient  to  sense  this  gradient  and  thus 
affect  the  accuracy  still  further.  A  full  quantitative 
evaluation  of  all  these  factors  and  how  they  Interact  and 
depend  on  probe  orientation  for  a  particular  flow  vrould 
require  estensive  tests  with  the  probe  immersed  in  flov/s 
whose  properties  \;ere  determined  in  some  alternate  manner. 
Such  tests  would  represent  a  sizeable  experimental  program 
in  themselves  and  so  none  have  been  attempted.  Application 
of  existing  boundary  correction  procedures  to  probe  data 
in  this  case  was  felt  to  be  of  little  value.  However,  it 
is  possible  to  draw  some  general  conclusions  about  the 
validity  of  the  measured  profiles. 

Tangential  velocity.  Effects  due  to  high  shear  should 
depend  on  whether  the  lateral  or  longitudinal  holes  are 
used  to  evaluate  this  component.  Check  runs  (figure  1?) 
using  both  methods  seem  to  bear  out  the  reasoning  that  the 


N  * -900 
L/D  =  5 
z/L  =  I 


V/U  =  1.5 


V/U  =  3.0 


20 


20h  <fi 


V,  15 


V, 


15 


o  X 


o  Axial  probe  travel 

(  lateral  hole*  sense  v) 

X  Transverse  probe  travel 
(  longitudinal  holes  sense  v  ) 


Figure  17.  Comparison  of  tangential  and  axial  velocity 
profiles  using  both  probe  traverse  methods. 


longitudinal  holes  nir.y  give  lower  values  than  the  lateral 
holes.  On  the  other  hand,  the  sol id -boundary  Interference 
effect  could  be  appreciable.  Even  at  a  distance  of  3 
4  probe  diameters  from  the  wall  (outer  40  per  cent  of  the 
profile),  such  an  effect  could  be  felt.  Again,  the  method 
of  measurement  should  be  a  factor.  It  appears  quite 
consistent,  therefore,  to  reason  that  the  outer  portions  of 
the  V  profiles  shovm  in  figure  15  exceed  the  Inviscld 
profile  due  to  v/all  effects  on  the  probe.  It  is  also 
consistent  to  note  that  the  check  runs  using  the  longitvidinal 
holes  to  determine  v  gave  lower  values  in  this  region. 

Axial  velocity.  A  higher  degree  of  precision  in 
measuring  axial  velocity  components  appears  to  be  inherent 
in  either  method  of  measurement  except  in  regions  where 
the  tangential  component  is  of  comparable  magnitude.  The 
close  agreement  shoT/n  by  the  check  runs  appears  to  bear 
this  out. 

Radial  velocity.  Although,  as  has  just  been  discussed, 
it  is  extremely  difficult  to  measure  tangential  and  axial 
velocities  vjith  any  great  degree  of  precision,  the  measure¬ 
ments  are  reasonably  reliable  quantitatively  and  are  certain¬ 
ly  accurate  enough  to  give  qualitatively  the  behavior  of 
these  profiles  with  variations  of  the  basic  parameters. 

On  the  other  hand,  since  the  radial  velocities  that  must  be 
measured  in  such  vortices  are  an  order  of  magnitude  smaller 
than  the  other  velocltjr  components,  the  measurements  that 
were  made  of  these  velocities  are  probably  of  very  little 
value  even  in  a  qualitative  sense.  Nevertheless,  for 
completeness,  as  \j'ell  as  to  give  some  idea  as  to  the 
magnitude  of  the  experimental  errors  that  were  Involved  in 
the  measurement  of  the  velocity  component  profiles,  a 
typical  measurement  of  radial  velocity  is  shown  in  figure 
18.  Included  in  this  figure  for  comparison  with  the  local 
measurements  are  the  radial  velocity  at  the  wall  as  obtained 
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Figure  18.  Typical  radial  velocity  measurements  compared 
with  the  profile  computed  from  the  measured  axial  velocity 
distribution. 
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by  mass  flov/  measurement  (the  solid  symbol)  and  a  curve 
showing  an  estimate  of  what  the  actual  character  of  the 
radial  velocity  profile  must  be  like.  This  es+^lmate  was 
made  by  assuming  that  the  axial  velocity  actually  satisfied 
the  equation  w  =  zf(r)  so  that  the  radial  velocity 
could  be  obtained  from  the  measured  axial  velocity  distribu¬ 
tion  through  Integration  of  the  continuity  equation.  It 
Is  seen  that,  as  might  be  expected,  the  order  of  magnitude 
of  the  radial  velocity  as  measured  Is  correct  but  the 
experimental  errors  are  such  that  no  Information  Is 
contained  In  these  data  concerning  the  character  of  the 
radial  velocity  profiles. 

6 A.  Static  pressure  profiles 

In  order  to  complete  the  picture  of  the  behavior  of 
flows  In  simple  cylindrical  vortex  chambers,  there  are 
presented  In  figures  19  and  20  data  concerning  the  static 
pressures  within  a  typical  vortex.  Figure  19  shows  typical 
plots  of  the  non-dlmenslonal  pressure  distributions  obtained 
across  the  exit  of  the  vortex  as  the  parameter  V/U  Is 
varied  for  the  specific  case  when  L/D  =  5  and  N  =  -900. 
Figure  20  shows  for  the  same  case  a  comparison  of  the  static 
pressure  at  z/L  =  0  and  at  z/L  =1.0  on  the  axis  of 
symmetry  as  a  function  of  V/U.  It  Is  Interesting  to  note 
the  change  In  the  character  of  the  static  pressure  distribu¬ 
tions  as  V/U  Is  increased  above  (V/U)  ,  which  in  this 

case  is  2.95^  and  as  V/U  is  Increased  beyond  a  second 
critical  value  at  V/U  -  I5  where  a  Jump  In  the  static 
pressure  on  the  centerline  takes  place  for  z/L  =  1. 


100 


Figure  19 .  Measured  static  pressure  distributions  across 
the  exit  plane  (z/L  =1). 


Figure  20.  Comparison  of  the  static  pressures  measured  on 
the  center  line  at  the  bottom  (S/L  =0)  and  top  (z/L  -  l) 
of  the  vortex  cup  as  V/U  is  varied. 
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7.  Discussion  of  Experimental  Results 

The  detailed  behavior  of  the  vortices  on  which  experi¬ 
mental  data  have  just  been  presented,  especially  the  way 
In  which  these  vortices  form  their  second  cell,  can  best 
be  understood  by  means  of  a  detailed  discussion  of  a 
typical  case.  We  v/lll  discuss  in  detail  the  vortex 
behavior  for  the  case  L/D  =  5  and  N  =  -900  with  the 
aid  of  figures  21  and  22.  Figure  21  consists  of  the  data 
already  plotted  in  figure  13>  together  with  the  follov^lng 
additional  Information:  (1)  theoretical  curves  for 
(p  -  2  versus  V/U  obtained  from  the  laminar 

theory  for  one -celled  (dash  dot  curves)  and  two-celled 
(dash  double  dot  curves)  vortex  configurations;  (2)  the 
pressure  drop  that  might  be  expected  if  the  fluid  at  the 
exit  of  the  vortex  were  turning  as  a  solid  body  (dashed 
curve);  and  (3)  the  pressure  drops  that  one  obtains  from 
an  integration  of  the  approximate  radial  momentum  equation 


dp  _  pv^ 
df  f 

using  the  measured  tangential  velocity  profiles  at  several 
values  of  V/U  (indicated  by  the  solid  symbols).  Figure 
22  gives  a  pictorial  representation  of  the  behavior  of  the 
streamlines  in  and  around  the  vortex  chamber  as  the 
rotational  speed  of  the  chamber  v;all  is  increased.  Also 
shown  in  this  figure  are  the  tangential  velocity  profiles, 
typical  of  each  condition,  in  the  exit  plane  of  the  vortex. 

Let  us  follov/  the  behavior  of  the  vortex  on  figures  21 
and  22  as  V/U  is  increased  from  zero  to  the  highest  value 
at  which  measurements  were  made.  Initially,  when  the 
tangential  velocity  is  zero,  the  flow  in  the  vortex  cap 
is  similar  to  that  shown  in  figure  22a.  The  axial  flow  in 


Ficure  Jl.  Compari3on  of  measured  pressure  behavior  at 
point  B  with  that  given  by  the  laminar  theory. 


Figure  ^^2.  Pictorial  representation  of  streamline  patterns  as  V/U  is  increased 
The  approximate  tangential  velocity  profiles  are  shown  in  the  exit  plane. 


Figure  22  -  concluded 
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the  exit  plane  is  everywhere  outward  with  a  velocity 
distribution  that  is  very  close  to  the  theoretical  curve 

p 

for  a  one-celled  laminar  vortex,  namely  w/w(0)  =  cos(f/R)  . 
There  is,  to  the  scale  of  the  static  pressure  plot  shown 
in  figure  21,  essentially  no  difference  between  the  static 
pressure  on  the  vortex  centerline  in  the  exit  plane  and  the 
ambient  pressure  (also  see  figure  19).  The  total  pressure 
shown  in  figure  22  for  V/U  =  0  is  higher  than  ambient 
pressure  by  Just  the  dynamic  pressure  of  the  outgoing  flow. 

As  V/U  is  Increased  from  zero,  the  flow  within  the 
vortex  cup  retains  its  one-celled  character  as  shown  in 
figure  22b  for  V/U  «  2.  Nevertheless,  the  dynamic  pressure 
of  the  exiting  fluid  on  the  centerline  of  the  vortex  is 
very  rapidly  reduced.  This  la  evident  both  from  figure  21 
and  from  the  actual  w/w(0)  velocity  profiles  for  this 
particular  case  which  are  shown  in  figure  15.  Apparently, 
a  second  cell  is  forming  somewhere  high  above  the  vortex 
chamber  and  is  gradually  moving  tovmrds  the  exit  plane  as 
V/U  is  Increased. 

We  should  also  note  in  figure  22b  that  the  entrained 
flow  is  no  longer  captured  in  Just  the  same  manner  as  it 
was  when  V  =  0.  This  is  due  to  the  centrlfvigal  pumping 
action  of  the  rotating  face  plate  of  the  vortex  chamber. 

As  a  result  of  this  pumping  action,  there  is  a  stagnation 
point  (actually  a  stagnation  ring)  somewhere  on  this 
rotating  plate.  Tlie  air  above  the  streamline  to  this 
stagnation  point  is  entrained  by  the  flow  out  of  the  vortex 
chamber.  The  air  below  this  streamline  is  forced  outward 
away  from  the  vortex  by  the  pumping  action  of  the  face 
plate . 

As  the  rotational  speed  is  Increased  further,  the 
dynamic  pressure  of  the  axial  flow  exiting  from  the 
vortex  becomes  less  and  less  until  at  the  critical  value 
(V/U)„_  =2.95  it  is  exactly  equal  to  zero  as  shown  in 

Cl* 
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figure  22c.  We  see  that  the  second  cell  that  had  developed 
far  from  the  vorte::  at  low  rotational  speeds  has  been 
sucked  down  toward  the  vortex  chamber  until,  at  V/U  =  2.95* 
the  tip  of  this  second  cell  has  just  reached  the  exit  plane 
of  the  chamber.  Up  until  this  critical  value  of  V/U, 
the  tangential  velocity  distribution  at  the  exit  plane  has 
not  been  greatly  altered  in  character  as  may  be  seen  from 
the  tangential  velocity  distributions  shovm  in  figure  15. 
Thus,  in  the  range  0  <  V/U  <  (V/U)  ,  the  pressure  on  the 

C/  X 

centerline  drops  uniformly,  roughly  as  the  square  of  V/U. 

Once  the  critical  value  of  V/U  has  been  reached,  the 
character  of  the  taiiLientlal  velocity  distributions  in  the 
exit  plane  changes  markedly  with  V/U  as  the  second  cell  of 
the  vortex  is  sucked  inside  the  chamber.  Indeed,  the 
tangential  velocities  are  lowered  so  much  in  the  center  of 
the  vortex  (see  figure  15)  that  the  pressure  drop  is 
actually  reduced  somewhat  as  seen  in  figure  21.  A  typical 
streamline  pattern  for  this  regime  is  shown  in  figure  22d. 

When  the  second  cell  of  the  vortex  has  been  sucked 
completely  into  tho  vortex  chamber,  v/hich  occurs  in  the 
case  under  discussion  at  V/U  =  8,  the  pressure  drop 
across  the  vortex  starts  to  increase  again.  It  will  be 
noted  however  that,  instead  of  follov/ing  a  curve  similar 
to  those  for  constant  Reynolds  nvimber  as  given  by  the 
laminar  theory,  the  pressure  drops  far  less  rapidly.  This 
behavior  in  all  cases  appears  to  be  due  to  an  effective 
increase  in  the  turbulent  or  RejTaolds  stresses  with  Increase 
in  V/U  as  discussed  in  Sections  3  and  4.  Streamlines 
typical  of  this  flow  regime  are  shown  in  figure  22e. 

After  the  two -celled  vortex  pattern  has  been  complete¬ 
ly  established,  further  Increase  in  V/U  results  in  a 
continued  Increase  in  the  pressure  drop  across  the  exit  of 
the  vortex  xintil  a  second  critical  value  of  V/U  is  reached. 
At  this  point,  V/U  ®  15  in  the  case  under  discussion. 
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there  is  a  rapid  rise  In  the  general  level  of  pressure 
within  the  vortex  at  the  exit  plane.  By  careful  probing 
of  the  flow  external  to  the  vortex  cup  with  fine  tufts 
and  with  smoke  streams,  it  was  possible  to  show  that,  in 
all  cases,  this  pressure  rise  was  associated  with  a  rather 
sudden  change  in  the  flow  pattern  outside  the  vortex 
chamber.  The  nature  of  this  change  is  shown  in  figure  22f. 
Inspection  of  the  streamline  patterns  for  V/U  <  15  in 
figure  22  reveals  that  there  is  alvjcys  an  entrained  flow. 
However,  as  V/U  Is  increased,  the  streamlines  leaving  the 
vortex  chamber  do  so  In  a  more  and  more  radial  manner. 

It  was  fovind  that  at  the  second  critical  value  of  V/U 
the  flow  leaving  the  vortex  cup  attached  Itself  almost 
Immediately  to  the  rotating  face  plate  outside  the  vortex 
thus  eliminating  the  entrained  flow  altogether.  The 
spreading  out  of  the  streamlines  outside  the  vortex  that 
resulted  was  found  to  greatly  reduce  the  amoiint  of  angular 
momentum  that  the  fluid  entering  the  second  cell  of  the 
vortex  could  obtain  by  viscous  action  from  the  primary 
flow  leaving  the  vortex.  As  a  result  of  this  lowering  of 
the  tangential  velocity  in  the  exit  plane,  the  pressure 
drop  across  the  vortex  was  greatly  reduced. 

It  would  be  expected  that  since  this  change  of 
streamline  configuration  must  depend  on  the  ratio  of  the 
general  level  of  axial  velocity  leaving  the  vortex  to  V, 
the  deeper  the  vortex  chamber  (higher  w)  the  higher  would 
be  the  V/U  required  to  effect  this  change.  The  fact  that 
this  is  so  and  that  all  the  vortex  chambers  investigated 
behaved  in  essentially  the  manner  described  above  is  amply 
demonstrated  by  figure  23  when  all  the  static  pressure 
measurements  given  in  figures  9  through  13  are  replotted. 

It  can  be  seen  from  figure  23  that  since  the  vortex 
cup  is  shallow  in  the  cases  of  small  L/D,  the  air  in  the 
second  vortex  cell  receives  far  less  angular  momentum  from 


Figure  23.  Summary  of  static  pressure  data  as  a  function  of 

V/U  and  L/d, 
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the  primary  vortex  flov/  than  In  cases  when  L/D  Is  large. 

This  is  particularly  true  for  the  case  when  L/D  =  1 
after  the  flow  outside  the  vortex  has  attached  to  the 
rotating  face  plate  (V/U  >8). 

Another  fact  that  is  demonstrated  by  figure  23  is 
that  V/U  Is  the  primary  parameter  governing  the  behavior 
of  the  vortices  under  discussion  here.  It  will  be  noticed 
that  for  eaoh  vortex  chamber  tested  (L/D  varied),  the 
roouits,  though  slightly  dependent  on  Reynolds  nimiber  are 
essentially  controlled  by  the  parameter  V/U,  This  Is 
particularly  true  of  the  vortex  chambers  of  low  L/D  when 
the  effects  of  molecvilar  viscosity  v;hich  are  expected  to 
be  largest  on  the  axial  flow  must  necessarily  be  small. 

The  results  shovm  in  figure  23  as  well  as  the  basic 
measurements  of  (V/U)  shown  in  figure  l4,  where  it  is 

C  I* 

seen  that  for  all  intents  and  purposes  (V/U)  is 
Independent  of  Reynolds  number  for  L/D  <  5#  give  strong 
support  to  the  basic  conjectures  of  Sections  3  and  4 
concerning  the  nature  of  turbulent  vortices. 

In  this  regard,  it  may  be  of  some  Interest  to  determine 
the  order  of  magnitude  of  the  effective  viscosity  that  a 
laminar  flow  (as  given  by  the  exact  solutions  of  Donaldson 
and  Sullivan)  would  have  to  have  in  order  to  match  the 
static  pressure  drops  that  have  been  observed  in  the  exit 
planes  of  the  several  vortex  configurations  tested. 

Figure  24  is  a  plot  of  the  ratio  of  effective  to  molecular 
viscosity  obtained  by  matching  the  experimental  static 
pressure  drops  at  the  exit  to  those  given  by  the  laminar 
theory.  This  matching  was  only  carried  out  when  fully 
developed  one-  or  tv/o-celled  vortex  structures  existed. 

The  data  are  not,  since  the  measured  flows  are  not  identical 
to  those  given  by  the  laminar  theory,  of  fundamental 
importance.  Nevertheless,  the  general  trend  of  the  data 


Figure  24.  Ratio  of  effective  to  molecular  viscosity  as  a 
function  of  V/U. 
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indicates  an  lncrer.se  in  effective  viscosity  which  Is 
roughly  proportional  to  V/U, 

One  final  point  should  be  made  before  bringing  this 
discussion  of  the  experimental  results  to  a  close.  This 
point  has  to  do  v/lth  the  self-consistency  and  accuracy  of 
the  data.  The  general  agreement  between  the  measured 
static  pressure  drops  and  the  pressure  drops  shown  In  figure 
21  computed  from  the  measured  tangential  velocity  distribu¬ 
tions  can  be  taken  as  a  measure  of  the  accuracy  of  this 
velocity  distribution  data.  Since,  in  general,  very  small 
changes  in  the  tangential  velocity  near  the  center  of  a 
vortex  cause  large  changes  in  the  static  pressure  drop 
across  the  vortex  it  may  be  concluded,  in  view  of  the 
general  agreement  of  the  results,  that  the  tangential 
velocity  distributions  measured  near  the  center  of  the 
vortices  in  this  study  are,  in  all  probability,  reasonably 
accurate . 
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8.  Conclusions 


As  a  result  of  the  experimental  study  of  the  flow  in 
simple  cylindrical  vortex  chambers  completely  open  at  one 
end  that  has  been  presented  in  the  preceding  sections,  the 
following  conclusions  may  be  dravm: 

(1)  For  turbulent  vortices  v/lthin  cylindrical  vortex 
chambers  having  length-to -diameter  ratios  In  the  range 

0  <  L/D  <  5»  the  primary  parameter  governing  the  character 
of  the  flow  is  V/U;  that  is,  the  ratio  of  a  characteristic 
tangential  to  a  characteristic  radial  velocity.  This  is  in 
marked  contrast  to  the  behavior  of  laminar  vortices  where 
in  addition  to  V/(J  the  radial  Reynolds  number  N  =  UR/v 
is  a  most  Important  parameter. 

(2)  For  turbulent  flow  in  vortex  chambers  in  the  above 
L/D  range,  the  establishment  of  a  two -celled  vortex  conflgura 
tlon  is  dependent  on  V/U  alone.  In  particular,  a  second 
cell  will  begin  to  be  established  in  any  such  vortex  when 
the  ratio  of  characteristic  tangential  to  radial  velocities 
exceeds  approximately  3* 

(3)  It  appears  that  turbulent  vortices  behave  in  an 
effectively  more  and  more  viscous  manner  the  higher  the 
ratio  V/U  becomes. 

(4)  The  experimental  results  tend  to  give  general 
support  to  the  conjectures  concerning  the  nature  of  three- 
dimensional  turbulent  vortices  given  in  Sections  3  and  4 
of  the  present  paper. 
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